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Nonlinear energy-conserving drift-fluid equations that are suitable to describe self- 
consistent finite-/? low-frequency electromagnetic (drift-Alfven) turbulent fluctua- 
tions in a nonuniform, anisotropic, magnetized plasma are derived from a variational 
principle. The variational principle is based on a drift-fluid Lagrangian that contains 
linear and nonlinear E x B velocities derived directly from the corresponding single- 
particle finite- gyrocenter Hamiltonian (in the zero-Larmor-radius limit). Covariant 
electromagnetic effects introduce a magnetic generalization to the standard ion po- 
larization density as well as introduce a new ion magnetization current, which are 
both missing from existing gyrofluid and drift-fluid Poisson- Ampere equations. An 
exact energy conservation law is also derived directly from the drift-fluid Lagrangian 
by application of the Noether method. 

PACS numbers: 52.30.Ex, 52.30.Gz 

I. INTRODUCTION 

The derivation of dissipationless nonlinear reduced-fluid equations that explicitly conserve 
energy represents an important goal of the numerical simulation communities in fusion [H-jst 
and space j^-Q] plasma physics. The existence of an exact energy conservation law plays 
an important role in benchmarking numerical investigations of energy-transfer processes in 
a nonlinear plasma turbulent state between the fluid kinetic energy, its internal energy, and 
the electromagnetic field energy. Hence, if and when dissipative effects and/or advanced 
fluid-closure schemes are introduced in a nonlinear dissipationless reduced-fluid model, one 
can develop a deeper understanding of the roles played by dissipationless and dissipative 
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effects on the long-time evolution of the turbulent plasma state. 

The variational formulation of dissipationless fluid dynamics ensures that the nonlinear 
fluid equations satisfy an exact energy conservation law derived by application of the Noether 
method Previous variational formulations of reduced fluid models include the drift- 



fluid work of Pfirseh and Conva-Rrstrcpo 1£ and the recent gyrofluid works of Strintzi 

12 1 - In contrast to Refs. Q,ll2|, however, we 



and Scott 1161 and Strintzi, Scott, and Brizard 



ignore finite-Larmor-radius (FLR) effects to focus our attention on the drift-fluid polarization 
and magnetization effects introduced by the low-frequency magnetic field fluctuations. 

In the energy- conserving drift-Alfven model presented here, the perturbed electric and 
magnetic fields 

E = -V0 -^1 and B ± = VAn X b , (1) 

c at 

are expressed in terms of the fluctuation scalar potentials <fi and A\\, and the time- 
independent background magnetic field B = B b is assumed to be spatially nonuniform. 
From Eq. (flj). we define the following perturbed E x B velocities: 

u E es |xV^ = E ± x^, (2) 



where Ug denotes the linear perturbed ExB velocity and V\\ denotes the parallel component 
of the nonlinear perturbed ExB velocity. Lastly, according to the finite-/? ordering [18]- 
[2^ (m e /mi < (3 1), we neglect the perpendicular component of the fluctuating vector 
potential (A^ = 0), so that the parallel component of the fluctuating magnetic field is 
assumed to vanish in what follows. 

The remainder of this paper is organized as follows. In Sec. |Hj we present the heuristic 
derivation of the nonlinear drift-fluid Lagrangian density to be used, in Sec. IIII1 in the 
variational derivation of our self-consistent, energy-conserving drift-Alfven equations. In 
Sec. IIV1 we use the Noether method to derive the local form of the energy conservation law 
for our drift-Alfven model. Next, we integrate the local form of the energy conservation 
law to obtain its global form, in which the energy-transfer processes between the parallel 
kinetic energy, internal energy, and electromagnetic field energy are manifestly expressed. 
In Sec. El we summarize our variational derivation of the energy- conserving drift-Alfven 
model. Lastly, in Appendix A, we present a procedure developed for including heat fluxes 
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within a fluid Lagrangian variational principle, while in Appendix B, we briefly outline the 
diamagnetic-cancellation procedure involving the addition of energy- conserving gyroviscous 
terms in the drift-fluid equations for (u\\,p±,p\\). 



II. DRIFT-FLUID LAGRANGIAN DENSITY 



The four-moment drift-fluid Lagrangian for the present anisotropic-temperature finite-/? 
electromagnetic model is based on the previous electrostatic gyrofluid work of Strintzi, Scott, 
and Brizard Here, the drift-fluid Lagrangian density (for a fluid species of mass m and 
charge e) is defined as 

L df = ^ u\ - (p ± + + en (A + A\\ b Q ) • " - en (f) - L ZLR [ H 2gy ] , (4) 

where u\\ = u ■ bo denotes the drift-fluid velocity parallel to the background magnetic field 
and all terms (except the last one) are standard fluid Lagrangian terms. The last term rep- 
resents the zero-Larmor-radius (ZLR) limit of the electromagnetic nonlinearities contained 

n 

in the second-order finite-/? gyrocenter Hamiltonian 21]: 

where v\\ denotes the parallel gyrocenter particle velocity, the E x B velocities ug and Vji are 
defined in Eqs. (J2J) and Q, respectively, and the background magnetic- field nonuniformity 
was neglected. By heuristically transforming this second-order gyrocenter Hamiltonian to its 
drift-fluid version, following a procedure outlined in Refs. the drift-fluid Lagrangian 

density © for each fluid species becomes 



mn 
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u\\ (b + B ± /5 ) +u E \ 2 - [p x + y) + en(A + A||b Q ) - en 



+ enA .» c - (p ± + ^) - en(<,* - U -iA;), (5) 
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where we introduced the definitions 



U\\ = Unb = m (1 + \B ± \ 2 /B 2 ) 1/2 ^ 



e$* = e<\> - m\u E \ 2 /2 



{e/c)A\ = (e/c)A\\ 



m V\\ 



(6) 



with the effective potentials $* and At both containing linear and nonlinear field terms. 
We point out that, while u\\ denotes the parallel drift-fluid velocity along the unperturbed 
(background) magnetic field B , U\\ denotes the parallel drift-fluid velocity along the total 
magnetic field B + B^. Moreover, although the term b 2 — 1 is considered small, it is kept in 
the drift-fluid Lagrangian (jSJ) to allow a covariant treatment of drift-fluid polarization and 
magnetization effects in our model [see Eqs. (fT0|) - |TT| below]. 

We note that the fluid velocity u appearing in the drift-fluid Lagrangian density (JSJ 



is defined as the particle-fluid velocity expressed in terms of moments o 
independent gyrocenter distribution function F(R, v», fi, t) defined in Ref. 



the gyroangle- 



n(r, t) u(r, t) = J /(x, v, t) v 5 3 (x - r) d 3 x d 3 v (7) 
= J F(R, v hf i, t) ((v ± + Vgc + • • •) 5 3 (R + p - r)) d 6 Z. (8) 

Here, (■ ■ •) denotes a gyroangle average and we used the transformation from a particle- 
fluid moment (JJJ) to a gyrofluid moment (JHJ), which was introduced in Ref. |23|, where v gc 
denotes the guiding-center (gc) velocity and higher-order corrections (e.g., electromagnetic 
fluctuations) are omitted for simplicity. Furthermore, because of the presence of 5 3 (R+p— r) 
in Eq. (JHJ), the gyrofluid moment of the perpendicular particle velocity vi does not vanish 
but, instead, yields the divergenceless term 



J F(R,v h n,t) (v ± 5 3 (R + p-r)) d 6 Z = - V X lp± 



cbo_' 



which introduces an important diamagnetic contribution to the particle flux 



cbn / „ , ^ r \ «_, / cb, 

/'Up : —prX (p± V In B + p|| b ■ Vb J - Vx lp L 



eB, 



o 



o 



eB, 



rb °x(V-P), (9) 



eB, 



o 



where P denotes the Chew-Goldberger-Low (CGL) pressure tensor 2J| [see Eq. (|21j) below]. 
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The appearance of the perturbed electric and magnetic fields in the drift-fluid Lagrangian 
(JHJ) yields the following expressions for the drift-fluid polarization and magnetization vectors: 

Px-^ = E^ F (Ex + Jb xB ± ), (10) 
9Ldf x - u\\ / Bj_ „ cb \ , \ 

M ^alI = 5:™^(«ii^ + Ex X ^j. (ii) 

Because of the mass dependence in the drift-fluid polarization and magnetization vectors 
([10 )1 - ([11 )1 . the drift-fluid polarization and magnetization effects are especially important for 
ion fluid species (in a quasineutral plasma). 

The total Lagrangian density for our four-moment electromagnetic drift-fluid model is, 
therefore, defined as the sum of the Lagrangian density of the electromagnetic field and the 
drift-fluid Lagrangian density (JHJ) for each fluid species (summation over species is implied): 

L = ±(|EJ*-|B,lV-( A °-" + 7^-**) 

+ ^l?-(- + f). < 12 > 
where we have omitted the contribution from the parallel electric field in the electric field 
energy (i.e., |E| 2 — > |E^| 2 ), which removes the parallel displacement current (dtEu) in the 
parallel drift-fluid Ampere equation [see Eq. (fHfiJ)]. The variational fields are the four drift- 
fluid moments (n, \i,p±,p\\) for each fluid species and the electromagnetic potentials (0, A\\). 
Here, we note that only the first four drift-fluid moments (n, u,p±,p\\) appear in the drift- 
fluid Lagrangian density © and higher-order drift-fluid moments (e.g., heat fluxes) are 
omitted as variational dynamical fields; future work will consider the introduction of higher- 
order drift-fluid moments in a Lagrangian variational setting j^j (see Appendix A for further 
comments). Lastly, we note that, because the background magnetic field Bo is independent 
of the variational fields, the magnetic-energy term — |B | 2 /87r has been removed from the 
Lagrangian density (fT2|) . 

III. VARIATIONAL DERIVATION OF DRIFT-FLUID DYNAMICAL 

EQUATIONS 

The drift- Alfven variational principle 

J 8h(n, u,p\\,p±; (f), A\\\ x) d 3 xdt = (13) 
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associated with the drift-fluid Lagrangian density (|T2|) must be expressed in terms of the 
Eulerian variations (Sn, Su, Sp±, Sp»; Scf>, SA^): 

( s dh 3h 3h . dL !± dh _ A <9L\ 

SL = [Sn — + Su- — + Sp\\ — + 5p ± — — + S<p — + SA\\ 



dn du dpn dp± dcj) dAn J 

where the additional x-dependence in L(- • • ; x) arises from the nonuniform background 
magnetic field B . Here, the Eulerian variations for the drift-fluid moments (Sn, Su, Sp±, 8p») 
are not independent of each other but are instead expressed in terms of a virtual fluid 
displacement £ (for each fluid species). 

The derivation of the Eulerian variations (Sn, Sp±, Sp\\) are based on dynamical constraints 
for the drift-fluid moments n a = (n,p±,p\\). The first dynamical constraint is associated with 
mass conservation, expressed in terms of the continuity equation 

y+fiV.u = fl, (15) 

where d/dt = d/dt + u- V denotes the total time derivative. The next two dynamical 
constraints are associated with the conservation of the first two single-particle adiabatic 
invariants, expressed in terms of the Chew-Goldberger-Low equations for the perpendicular 



and parallel pressures 



^ + p ± V-u + p± (I — b b ) : Vu = 0, (16) 
-j£ + p|| V-u + 2p|, b b : Vu = 0, (17) 

where the higher-order heat-flux moments and are omitted here but are considered 
later in Appendix B. 

We begin our derivation of expressions for the Eulerian variations to be used in Eq. ()14j) 
by transforming the dynamical constraints (|15|) and (|16 )) -()17 )) into Lagrangian variations 
Ar] a = (An, Ap±, Apn) by introducing the following limiting process 

lim | At — | = An a and lim (uAt) = £ (18) 
At-*o \ dt ) At^o v ' v ' 

where ^ denotes the virtual fluid-displacement (for each fluid species), with the Lagrangian 
variation of the fluid velocity is defined as Au = d^/dt. Hence, using the procedure (|18j) on 
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the constraint equations (|15 |1 -(|17 |1 . we find the Lagrangian variations 



An = -jiV-d 
Ap ± = -p±V-£-p± (I-b b ):V£, 
Ap|| = -pi, V-£-2p ± b b : V£. 

Next, using the relation 5r] a = An a — £ • Vr] a between Lagrangian and Eulerian variations 
2fi| . the Eulerian variations (5n, 5u, Sp\\, 5pj_) to be used in the variational principle (|T3j) are 
now defined as 

Sn = - V • (n £) 



5u = d£ + (u-V)£ - (£-V)u 



S P]l 



£) - 2pn b b : V£ 



(19) 



<fy± = - V • (p ± £) - P± (I - b b ) : V£ _ 
The remaining drift-fluid equations include the evolution equation for the parallel drift- 
fluid velocity u\\ and an explicit expression for the perpendicular component of the drift-fluid 
velocity u. Lastly, a self-consistent treatment involves the appropriate Poisson and parallel 
Ampere equations for the potentials (<j),A») expressed in terms of the drift-fluid moments 
(n,u\\). 



A. Variation of the Lagrangian density 



We now insert the Eulerian variations (fiT)|) into the drift-fluid Lagrangian variation ()14j) . 
and rearranging several terms to isolate £ and Sip 1 = (5<ft, 5 An), to obtain 



5L 



+ 6ip 



d (dh\ „ ( dL\ 



Vu 



dh 

du 



V-P 



rf V 



dr] a 



dh 

dip 1 



- V 



dh 



d(y*i) 1 



+ « +v - r ' 



where 



a „dh Bh flL flL 

V V 17T = ^ V- — hp,, V- hpx V- — 

c/r7 a on ap\\ op± 



(20) 
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and the tensor P denotes the CGL pressure tensor 



24] 



P = - P± 7^- (I - b b ) - 2p\\ — b b = p± (I - b b ) + p\\ b b . 



(21) 



While the space-time divergence terms <9fA + V ■ T in Eq. (1201) . expressed in terms of the 
Noether fields 



A 



and 



<9L 



<9(V# 



du 



+ 



i he 



(22) 



(23) 



do not play a role in the variational principle / SL d 3 x dt = 0, they p 
derivation of exact conservation laws based on the Noether method 



av a crucial role in the 



111 



'see Sec. Ha. 



B. Drift-fluid velocity equations 

The stationarity of the action functional / L d 3 x dt with respect to a arbitrary virtual 
fluid displacement £ yields the Euler-Poincare equation j3] (for each fluid species) 



d fdh\ „ ( dV 



Vu~ + V-P 
ou 



rf V 



(24) 



which describes the time evolution of the drift-fluid velocity u. Upon substituting derivatives 
of the drift-fluid Lagrangian density (fT2|) into the Euler-Poincare equation (J24|) . and using the 
fact that the background magnetic field B is assumed to be a time-independent nonuniform 
vector field, Eq. (J23|) can be written as 



= n b — mm b 2 + - At) uxB' + V-P + tiV 

at \ c V c 

where we have introduced the following divergenceless field 



m 



e$* + y(M||&) 2 



(25) 



B* = B + V x 



A * 7710 7 5 

At H u\\V 

e " 



(26) 



from which we define 



B\ = b -B* = B + 
,* B* , cb 



A\ + ) ,„ /r 



b • V x b , 



BT, 



eBl 



X 



mu\\ b 2 + -At)K, - V \mu\\b 2 + - A] 



(27) 
(28) 
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where k, = b ■ Vbo denotes the background magnetic curvature. We note that the term 
b • V X b in Eq. (|2~7|. which is related to the unperturbed parallel current flowing along the 
background magnetic field, may vanish for some magnetic geometries (e.g., magnetic-dipole 
geometry) . 

The Euler-Poincare equation (|2*5j) can be divided into two equations: one equation that 
expresses the drift-fluid velocity u in terms of the drift-fluid moments (n, u\\,p±,pn) and the 
scalar potentials ((f), An), and one equation that describes the time evolution of the parallel 
drift-fluid velocity u\\. 



1. Drift- fluid velocity 



The first equation is obtained by taking the cross-product of Eq. (f23j) with bo, which 
yields the following first-order expression for the drift-fluid velocity: 



" = «||b* + -^X 
enBu 



m 



V-P 4 // V (( •!• + y (u\\b) 2 



(29) 



where b* is defined in Eq. (|2*8)l. Here, the divergence of the CGL pressure tensor (|2"Tj) is 



V • P = V P± + p A + [V • (p A b ) 



(30) 



(with pa = p\\ — p± denotes the pressure anisotropy) so that b X (V • P) = b X (Vp± + 
Pa^o)- According to Eq. (j2HJ), the drift-fluid velocity u, therefore, consists of the paral- 
lel velocity u\\ = u • bo (note that bo • b* = 1) and the following perpendicular drift-fluid 
velocities 

u D = (cbo/enBjj xVpi 



u$ = (cb /eB^j x V (e$* +mC/f/2 



(31) 



u c = (c bo/en B*A xj?a «o 
corresponding to the diamagnetic velocity, the generalized E x B velocity (which includes 
the gradients of | | 2 and £7| = u"jb 2 ), and the curvature-drift velocity, respectively. We 
note that, although the higher-order polarization drift velocity 



u 



pol 



bo_ duE 

Q dt 



(32) 
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is not included in the drift-fluid velocity (j2HJ), it is self-consistently introduced in the net cur- 
rent density J (through the polarization current density c^PjJ in the form of the quasineu- 
trality condition V ■ J = (see Sec. IIIIDI for more details). 



2. Evolution equation for u\\ 



The second equation to be extracted from the Euler-Poincare equation ([25)1 . which de- 
scribes the time evolution equation for parallel drift-fluid velocity u\\, can be obtained by 
taking the dot-product of Eq. (|23jl with b*: 



9 f ,2 e A* 

n — mu\\ o + - A\ 
at \ 11 c 



-b*- 



m 



V-P + // V (< <!>' - '—XJ\ 



(33) 



Here, the terms on the left side are 



, du \\ 

+ Mil 



dt 



'V±A|| 



dA\\ 

~dt 



and 



e dA\\ 



— m 



dV\\ 
~dt 



c dt c dt 

which contains both the parallel inductive electric field and the partial time derivative of the 
parallel component of the nonlinear E x B velocity (jHJ); note that the parallel component 
of the partial time derivative of the linear E x B velocity vanishes (bo • dt\iE = 0). Note 
also that, although the convective derivative u ■ V(mw|| b 2 ) along the drift-fluid velocity u 
appears to be absent in Eq. ()33|). it is actually hidden in the last term —mnU\\ b* • VC/ii; 
this term may appear explictly in Eq. (|23jl by redefining the divergenceless vector (|2*Bj) into 
a non-divergenceless vector as was done in Ref . ^ | . 

The set of drift-fluid equations of motion for the four drift-fluid moments (n, u\\,pj_, p\\) 
are, thus, given by Eqs. (|I5 )) .(|33 |L and (|T^j) - (fT7j) . respectively, while the drift-fluid velocity u 
is defined by Eq. ()29j) . Note that each of the drift-fluid equations for (u\\,p±,p\\) involves the 
diamagnetic advection operator ud • V, which must be cancelled by the addition of so-called 
gyroviscous cancellations (see Appendix B for further details). Furthermore, the drift-fluid 
moment equations now need to be complemented by suitable evolution equations for the 
electromagnetic scalar potentials (<f>,An) through the low-frequency version of Maxwell's 
equations, in which the net charge and current densities are expressed in terms of drift-fluid 
charge and current densities as well as polarization and magnetization terms. 
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C. Poisson and parallel Ampere equations 

We now derive the drift-fluid Poisson and parallel Ampere equations from the variational 
principle (JTHJ) as Euler-Lagrange equations for ip l = ((j),A\\): 

dh dh 

= dfi ~ v " 9(v^)' ^ 

where the drift-fluid polarization and magnetization effects, represented by the vectors (jlOj) - 
flTTj) . are introduced naturally in these respective equations. 

The drift-fluid Poisson equation for <fi can be found from the Euler-Lagrange equation 
(J34j) for ip l = (p, which can be expressed as 



V-Ej 

An 



en + V± 



mnc 2 /_ , mii 



= E«»" V-P ± , (35) 

where summation over fluid species is shown explicitly. This equation shows the effects of 
the finite-/? electromagnetic generalization of the polarization density p po \ = — V • Pj_, where 
the drift-fluid polarization vector is defined in Eq. pOj) . 

Next, the parallel drift-fluid Ampere equation can be found from the Euler-Lagrange 
equation f° r ^ — which can be expressed as 
V*A\\ ^ en 

V c 



47T C 



^ en _ 



= E + V- (M ± xb ), (36) 

where the absence of the parallel displacement current dtE\\ results from the omission of in 
the drift-fluid Lagrangian density ()12j) . We also see that the second term in the divergence 
on the right side of the drift-fluid Ampere equation represents the effects of the parallel 
component of the magnetization current bo ■ J m ag = c b • V X Mj_, where the magnetization 
vector Mi is defined in Eq. (II 1J) and background magnetic nonuniformity is ignored. Since 
the polarization vector (jlOj) is perpendicular to the background magnetic field, however, the 
drift-fluid polarization current J po i = dPj_/dt has no parallel component (bo • J po i = 0) and, 
thus, it does not appear in the parallel drift-fluid Ampere equation (JHEJ)- 

Lastly, we note that the covariant electromagnetic polarization and magnetization effects 
in the drift-fluid Poisson- Ampere equations are due to (a) the — A11 symmetry of the nonlin- 
ear Lagrangian term mnun Vn in the drift-fluid Lagrangian (|5|L and (b) the nonlinear kinetic 
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energy terms mn(|u£;| 2 + L/|)/2. While such low-frequency polarization and magnetization 
effects appear naturally in nonlinear finite- (3 gyrokinetic theory I21I, th ey have not appeared 
explicitly in previous derivations of gyrofluid equations (e.g., Ref. |23|). 



D. Quasineutral drift-fluid dynamics 

A fundamental aspect of drift-fluid dynamics in a strongly magnetized plasma involves 
the fact that the time evolution of the turbulent plasma state maintains a condition of 
quasineutrality [2fl| . where p = V • Ej_/47r ee and, thus, according to the drift-fluid Poisson 
equation (f33|) . the net drift-fluid charge density [14| 

en = V • P ± (37) 

is expressed in terms of the drift-fluid polarization vector (fTUj) . This drift-fluid quasineutral- 
ity condition, therefore, can be used with Eq. (|33j) to obtain the drift-fluid charge conserva- 
tion law: 

= V- (j^enu + = V-J, (38) 

where the net current density 

jE E efiU + y + Vx X (39) 

is defined as the sum of the net drift-fluid current density (J2 en u) , the drift-fluid polariza- 
tion current (d t P±), and a divergenceless term (V X x)- 

The vector field x i R Eq. (|3*^j) is determined by considering the parallel drift-fluid Ampere 
equation (J36j) and introducing the definition for the net parallel current density: 

J,, ee -iLVjAii = V en Mil + cbo-VxMi, (40) 

4:71 

where the drift-fluid magnetization vector Mi is defined in Eq. (jllj) . Hence, we find x = 
c Mj_ and the net current density is now expressed as [l^ 

dP 

J = Y en u H — - + cVxMi 

at 

^ 

= J|]b + J2 enu ± + ~Qf + c(VxM ± ) ± , (41) 

where we have introduced the net parallel current density (}4T)j) and (V X Mj_)j_ denotes 
perpendicular components of the drift-fluid magnetization current. 
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Lastly, we note that, while the standard polarization drift velocity ()32|) is absent from 
the drift-fluid velocity (j2HJ), the associated drift-fluid polarization current density d t P± is 
introduced naturally into our formalism through the drift-fluid quasineutrality condition 
(jHZj) and the drift-fluid charge conservation law ()38j) . Furthermore, we note that the drift- 
fluid polarization vector (fTT)|) can be expressed as Pj_ = enw± so that the drift-fluid 
polarization current density 



dt ~ ^ 



en — : — — V • enu w i 

dt v ; 



can also be expressed in terms of the total time derivative d/dt = d/dt + u- V and the 
polarization velocity (}3*2"j) now appears explicitly in the term dw±/dt. 

IV. LOCAL AND GLOBAL DRIFT-FLUID ENERGY CONSERVATION LAWS 

In this Section, we present the local and global forms of the energy conservation law, as 
they arise from an application of the Noether method For this purpose, we point 

out that, as a result of the variational principle / SL d 3 x dt = 0, the only remaining terms in 
the Eulerian variation of the Lagrangian density (|2()jl yield the Noether equation 

5L = ^ + V-r, (42) 

where the Noether fields A and T are defined in Eqs. (|22j) - (J23|) . The energy and momentum 
conservation laws are derived from the Noether equation ()42j) by considering infinitesimal 
time and space translations Q|, respectively. In the present work, we focus our attention 
on the local and global energy conservation laws associated with our drift- Alfven model. 

A. Local energy conservation law 

We derive the local form of the energy conservation law from the Noether equation (}4*2"j) 
by considering infinitesimal time translations t — > t + St, from which we obtain the following 
expressions for the virtual fluid displacement £ and the Eulerian variations S<p and <5L: 

i = -uSt and (50, 6A h SL) = -St ^ . (43) 
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Inserting these expressions into Eqs. 



the Noether fields become 



A 3 -«(„•£) 



r = - st 



dL 



dL 



U U " 7T~ — V ~n 



+ P-u + 



dip* dL 



(44) 
(45) 



dt d(Vip' 

where summation over fluid species is, henceforth, implied wherever appropriate. By com 
bining Eqs. (|43|) - (|45j) . we obtain the primitive form of the local energy conservation law: 

de' 



dt 



+ V-S' = 0, 



(46) 



where the primitive energy density is 
dL 



£ = U 



du 



L 



2 11 V 2 
and the primitive energy-density flux is 



-nmU\\ + (p_L + ) + en$* + — (|V±Ay 



(47) 



S' 



<9L 



<9L 



P-u + 



dip 1 dL 



Pll 



9* d(V# 



u + P-u + 



(48) 



dt <9(V#) 

Note that the energy conservation law (|%5Jl has the following gauge-invariance property: the 
energy conservation law is unaffected by the transformation 



e = e' + V-D and S = S' 



<9D 

~dt : 



(49) 



where D is an arbitrary vector field. 

In order to arrive at the final form of the local energy conservation law, we need to 
rearrange some terms in Eq. (|47j) . First, by substituting the drift-fluid Poisson equation 
(J35J), we find 

|V0| 2 |V0| 2 mn 



en 



in 



<77 + 2 



\ u e\ 



mnuw V\\ 



VD, 



(50) 



where the gauge vector field D is defined as 

, dL 



D 



<9(V0) 



(51) 



15 



Next, we express the primitive energy density (|47|) as e' = e — V • D, where the final form 
of the energy density is defined as 



mn 



uii (b + BJB ) + u E 



+ (P± + 



+ IB 



at + v ' s 



The final form of the local energy conservation law is, therefore, expressed as 

0. 

where the final form of the energy density flux, defined as 

S 



q, dB 

s " IT' 



is expressed, after some partial cancelations, as 

'm 



n 



P± 



2 



u + P-u 



where the field energy density flux S„ is defined as 







dAu dL 



(52) 



(53) 



(54) 



(55) 



(56) 



dt \d(V(f)) J dt d(VA\\)' 
In Appendix B, after gyroviscous diamagnetic cancellations are inserted into the drift- 
fluid equation for (un,p\\,px), the local energy conservation law (jH^j) is converted into a new 
energy equation 



dt 



+ V-S* 



0. 



(57) 



in which diamagnetic-cancellation terms result in a heat-modified energy density flux S* [see 
Eq. IjBSjl ]. This new form ensures that the total energy £ = J ed 3 x still satisfies the global 
energy conservation law dS/dt = 0. 



B. Global energy conservation law 

The global energy conservation law 

can be derived from the local energy conservation law by integrating it over space (and 
neglecting surface terms), where the component energies 

£ u = J d 3 x t/f) , (59) 
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£ 

p 

Sep 

S A 



d 3 x (p ± + |) 
/ |V0| 2 



d 3 x 
d 3 x 



mn 



8tt 
|V±A|, 



+ — - |u E | - mnu\\ V\\ 



in 



(60) 
(61) 

(62) 



represent the parallel kinetic energy, the internal energy, and the electromagnetic field en- 
ergy (defined as the sum of the 0-field energy and the Ay-field energy Sa), respectively 
We now present the time evolution of each energy component (in the absence of diamag- 
netic cancellations and heat fluxes) in order to identify the energy-exchange processes that 
allow the transfer of energy between the three types of drift-fluid (parallel kinetic, internal, 
electromagnetic field) energies. 

First, we begin with the time derivative of the parallel kinetic energy density 







" (mn 2 \ 




= - V- 









( 1 dh 

nu\\ mu\\ o — + 



- u- (V- P + enV$*) 
e dAt\ 



(63) 



dt c dt t 

derived from Eq. ()33|) after rearranging some terms and using the continuity equation (JT5j) 
Hence, the time evolution of the parallel-kinetic energy is 



dt 



d x 



( db e dAl 
u • (V • P + enV$*) + nu\\ \ mu\\ b — H 



dt 



c dt 



(64) 



Next, we evaluate the time derivative of the internal energy density V = p± +p\\/2 using 
the parallel and perpepdicular CGL pressure equations (fI5 jl -(|I7 jl : 



dV 

-p + PV-u+P:Vu 

at 



0. 



so that, by rearranging terms, we obtain 

- V- (Pu + P-u) + u- (V-P) 



dV 
~dt 



and the time evolution of the internal energy is 



d£ p 
~dt 



J d 3 x u • (V • P) . 



(65) 



(66) 



Thirdly, we evaluate the time derivative of the Ay-field energy density 



d_ ( \V±A 
dt \ 



1 dAw 



71 



4tt dt 



VI An + V' 



'dAw ViAn 



dt An 
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- V 



dA\\ dh 



dt dVA 



en dA\\ 

+ u \\ 

c at 



mnu\\c 
R 2 



U\\ 



V id) V i A 



dA \\ 
~dt 



(67) 



where the second equality is obtained by substituting the drift-fluid Ampere equation (J3BJ), 
so that the time evolution of the An-field energy is 



dS A 
dt 



d 3 x 



( db e dA\ 

nu\\ mu\\ b — + — - 

11 \ 11 dt c dt 



mnu\\c 

n 



Lastly, we evaluate time derivative of the 0-field energy density 



0_ 

dt 



mn 







ITT 



+ We\ — mnu\\ V\\ ) = — en$* 

2 at \ 



57T 



V- - uen$* 

V dt 



+ en 



d<5>* 
~dt 



dA \\ 

In 



+ V-Dj 



(68) 



(69) 



where the first equality is obtained from the definition (J5Uj) while the second equality arises 
from the use of the continuity equation (|15jl. By using the definition 



en 



d$* 
~df 



V±0 y d<j> 
4tt ± dt 



dt 
- v 



en + V 



V(f> mnc 2 „ 

1^ + ^T V ^ 



-is 



dt V 4tt 



1st 



mnc 

R2 
■°0 



mnu\\ 



so that Eq. (|69|) becomes 
IV,^| 2 



at 



+ 



mn 



+ U- (enV$*) — mnw 



mnuw V\\ 



V i A 



d_ ( dh \ 
dt [ lN4> I 



en $* u 



(70) 



and the time evolution of the 0-field energy is 
dE<h 



dt 



d A x 



mnuw 



v i A,, 



(71) 



By combining Eqs. (|HHj) and (fTTj) . and using the expression 



dV\\ 
~dt 



R2 
-°0 
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with the definition © for A», we obtain the time evolution of the electromagnetic energy 



dt 



d 6 x 



( db e dA\ 
enu • V$* + nu\\\mu\\b— + 



dt 



c dt 



(72) 



We now see that by combining Eqs. (|b4|). (JBBj) . and (|72jl . the global energy conservation law 
(J55J) is satisfied exactly. Hence, the energy-transfer processes (jMJ), (JUSj), and (J72J) between 
the parallel kinetic energy £ u , internal energy £ p , and electromagnetic field energy 8^ involve 
the following pathways: 



(P±,P\\) ^ ( 73 ) 

where Eqs. (I64|) and (|66p. for example, show that internal energy S p is transferred to and 
from parallel kinetic energy S u . 



V. SUMMARY 

The nonlinear drift- Alfven equations (|15p. (jHHj) . and ()16| l -()17 p for the anisotropic drift- 
fluid moments (n,u\\,p±,p\\), respectively, and the drift-fluid Poisson- Ampere equations (}35j) 
and (jHfjjl were derived in this paper from a Lagrangian variational principle based on the drift- 
fluid Lagrangian density (|12|1. The local and global forms of the exact energy conservation 
law for these drift- Alfven equations, given by Eqs. and respectively, were derived 
by the Noether method as a natural consequence of the variational formulation presented 
here. 

The inclusion of the linear and nonlinear E x B velocities ®-(jHJ) as well as the parallel 
flow along the total magnetic field lines into the drift-fluid Lagrangian density (fT2*|) enabled a 
derivation of generalized (electromagnetic) drift-fluid polarization and magnetization effects. 
Under the drift-fluid quasineutrality condition, the drift-fluid Poisson- Ampere equations 
and (jHUj) were transformed into the drift-fluid charge conservation law (|38|). with the net 
current density J defined by Eq. (|4ip. which represents a generalized (electromagnetic) form 
of the so-called vorticity equation. 

Although higher-order drift-fluid moments have been omitted from the present variational 
formulation of our four-moment nonlinear anisotropic drift-Alfven model, Appendix A shows 
how higher-order heat fluxes might be introduced into a variational formulation, while Ap- 
pendix B shows how diamagnetic cancellations in the drift-fluid equations for (u\\,p±, pu) 
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can be incorporated without affected the global energy conservation law. Future work will 
investigate the extension of variational principles by including higher-order heat fluxes as 
dynamical variational fields as well as the extension of Generalized Lagrangian Mean (GLM) 
methods (^J to plasma fluid models. 
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APPENDIX A: NOTES ON MACROSCOPIC PLASMA LAGRANGIANS 

In this Appendix, we present an adapted version of Y.-K. M. Peng's Ph. D. dissertation 
j^j] (with selected parts also published in Ref. [ssj]). The macroscopic Lagrangian presented 
by Peng is based on constrained variations of the fluid density, velocity, and pressure tensor, 
and includes the effects of finite heat fluxes self-consistently. 



1. Lagrangian Density in Eulerian Coordinates 

Following Peng's work, we begin by combining the pressure constraint equations pfi|1 - (|17j) 
to form the scalar pressure equation for V = p± + p\\/2: 

dV 

— + PV-u+P:Vu = -V-q, (Al) 
at 

where the heat flux q now appears on the right side of Eq. (|A1|) . Next, the heat flux vector 
q is represented as 

where Q denotes the total density of heat energy transported across the boundary of the 
macroscopic (coarse-graining) cell as it moves along its Lagrangian trajectory 0,1^. Hence, 
the right side of Eq. (|A1J) can be written as 

dQ dV-Q 
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so that Eq. (|A1|) may be rewritten as 

^-{V + V-Q) = - PV-u - (P + VQ) : Vu. (A3) 

By using the Lagrangian- variation limitting procedure (j!8|) . we, thus, obtain the Lagrangian 
variation 

A(? + V-Q) = - PV-C- (P + VQ) : V£, (A4) 
from which we obtain the Eulerian variation 

6{V + V-Q) = - V-{P£ + £ • VQ) - P : V£. (A5) 

Lastly, the heat flux q is introduced in the Lagrangian density L, defined by Eq. (J12)) . by 
adding the divergence term —V-Q: 

L' = L - V-Q = L cold - (V + V • Q) , (A6) 

where \j co u denotes the cold-fluid part of the drift-fluid Lagrangian density ((T2J) • 

2. Eulerian Variational Principle 

We now write the expression for the Eulerian variation of the heat-modified Lagrangian 
density ifATSjt : 

5L' = 5L cold + V • (V £ + £ • VQ) + P : V£ = 5L + V • • VQ) , (A7) 



where the Eulerian variation SL is given by the Noether equation (J42J). We note that the 
quantity Q does not play a dynamic role since it appears inside a spatial divergence. We 
shall see, however, that it plays an important role in terms of the energy conservation law 
derived from the heat-modified Noether equation 

51! = ^ + V- (r + £-VQ). (A8) 

3. Energy Conservation Law by Noether Method 

The invariance of the Lagrangian density (|A6|) under time translations t — > t+St yields an 
energy equation derived from the Noether equation (|A8[) . By applying the Noether method 
(as outlined in Sec. llVjl . we obtain the heat-modified energy equation: 

= | (e + V-Q) + V- (S + u.VQ) s | + V (S + f ), (A9) 
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where the energy density e and energy-density flux S are defined in Eqs. (J52J) and (J3HJ), 
respectively. Lastly, using the definition (jA2|) for the heat flux q, we obtain the energy 
equation 

^ + V-S = - V-q. (A10) 

Although the heat flux q appears self-consistently within Peng's variational principle, we 
note that the evolution equation for q is outside of its scope. It is a topic of future research 
to determine whether the evolution equation for the heat flux can be derived within a 
variational formulation. 



APPENDIX B: DIAMAGNETIC CANCELATIONS AND ENERGY 

CONSERVATION 

The drift-fluid equations (|33|) and (|16|) - (jl7|) contain diamagnetic advective derivatives 
Ufl ■ Vi] a of the drift-fluid moments rf = (u\\,p±,p\\) that must be cancelled by the addition 
of so-called higher-order gyroviscous cancellations. Since these higher-order moments cannot 
be derived from the drift-fluid Lagrangian, the diamagnetic cancelations must be inserted 
in drift-fluid equations in the post-variational phase. However, there exists a constraint in 
the addition of higher-order moment terms, namely, that the global energy conservation law 
should not be altered. 



1. Parallel drift-fluid dynamics 

The diamagnetic cancelation needed for the parallel momentum equation ()33|) involves 
the addition of the term —V ■ II| on the right side of Eq. ()33|) . associated with the non- 



diagonal part of t 
is found to be 133 



re pressure tensor. From Vlasov theory, the diamagnetic-cancelation term 



lc press 



- V-IIf = -V- (p±^?xVm||) = mnxi D -Vu\\ + p ± K(mu\\) , (Bl) 

ei>|i 

where the term Bt replaces the denominator Bo to ensure an exact diamagnetic cancellation 
and the magnetic differential operator JC(- ■ •) is defined by the identity 

J^xV/ = ^-V/xVj - g)C(f), (B2) 



22 



valid for arbitrary functions / and g. 

Note that the diamagnetic-cancellation term (jBl|) is energy conserving since 

- u|| v-n; = - v - (u|| n|) 

and, thus, the time evolution of the parallel kinetic energy density (J53j) is now expressed as 

!(in If) =-V.(=pt7f« + «,ni) -u.(V.P + enV*') 

db e dAi\ 

+ (B3) 

and the global time evolution equation (JMj) for £ u is still valid. 




2. Internal energy 

To consider the diamagnetic cancellations in the anisotropic pressure equations ([16 p -([17j h 



we write the modified CGL pressure equations 



llp± + p ± V • u + P± (I - b b ) : Vu = - V • - 2 • (b • Vb ), (B4) 



Jl + p{] V .u + 2p n b b : Vu = -2V-q? + 4 q? } • (b • Vb ), (B5) 



dt 

dt 

which have been modified by the addition of the parallel and perpendicular gyroviscous heat 
fluxes q^ and q^ 1 , respectively, derived from the Vlasov equation directly and found to be 

(||) 1 cp±b (±) cp ± b 

q± = g ~^B*~ X VT " + p W Uc and q ^ 2 X VT± ' ( B6 ) 



Inserting these diamagnetic heat fluxes into Eqs. (jB4|) - (jB5|) . the time evolution of the internal 
energy can, therefore, be expressed as 

dV 

— = - V- (Vu + P-u + q) + u- (V-P), (B7) 
at 

where the net heat flux is q = q^ + q^. Hence, we see that the insertion of diamagnetic 
cancellations into our finite-/? electromagnetic drift-fluid model does not jeopardize the local 
and global energy conservation laws, where the energy density flux (J55|) is now replaced by 
the flux 

S* = S + uii n: + q? + q^. (B8) 
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